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Abstract 

A self-consistent Thomas-Fermi-Poisson based calculation scheme is used to achieve 
spin resolved incompressible strips (ISs). The effect of exchange and correlation is 
incorporated by an empirically induced g factor. A local version of the Ohm's law de- 
scribes the imposed fixed current, where the discrepancies of this model are resolved 
by a relevant spatial averaging process. The longitudinal resistance is obtained as 
a function of the perpendicular (strong) magnetic field at filling factor one and 
two plateaus. Interrelation between the ISs and the longitudinal zeros is explicitly 
shown. 
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Devising spin degree of freedom as a coherent state to carry information be- 
came an essential phenomena in today's quantum information processing re- 
search, also in the integer quantized Hall effect [1] (IQHE) regime. However, a 
consistent picture describing the current distribution within these samples is 
still under debate. Recent attempts to calculate these quantities, can be cat- 
egorized under two essential titles: (i) A semi-classical solution of the Poisson 
and Schrodinger equations to obtain electrostatic quantities [2,3], i.e. a Hartree 
type (self-consistent) Thomas- Fermi(-Poisson) approximation (SCTFPA), in 
which the imposed current is treated within a local version of the Ohm's 
law [4,5] (LOL), (ii) A full quantum mechanical solution of the Hartree- Fock 
Hamiltonian, either using direct diagonalization[6] or local spin density ap- 
proximation (LSDA) plus density functional theory (DFT) [7] . The latter two 
approximations, although much powerful to describe the electrostatics, lack 
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Fig. 1. The local filling factors (solid lines) calculated at two representative magnetic 
field values, ft = fttu c /Eo, where Huj c = HeB/m is the cyclotron energy and Eq 
is the pinch-off energy given by 2-7re 2 no/K, with the fixed donor number density 
(no = 4x 10 11 cm -2 ). Corresponding current densities (broken lines) a small current 
(7 = Uh/O, = 0.01, see Ref.[4] for the definition) is driven in y direction, which is still 
in the linear response regime. The calculations are done at the default temperature 
ksT/EQ = 0.001 for a sample width 2d = 6.3 /im. 

presenting the global longitudinal or Hall resistances, since the direct diagonal- 
ization model is essentially focused on describing the local compressibility or 
conductivity and the second one uses the commonly used Landauer-Biittiker 
type conductance phenomena. On the other hand, the simpler approach, i.e. 
SCTFPA+LOL, was able to provide information of the local current distri- 
bution and the global resistances as a function of magnetic field considering 
spinless electrons [8]. Here, we extend our previous works [5,8] to include the 
effect of spin, i.e. Zeeman splitting, by considering an experimentally induced 
g factor [9]. The simplest, but not the most trivial, treatment of Zeeman split- 
ting, conductivity model and exchange-correlation effects is used to describe 
the mentioned quantities. We, essentially incorporate Zeeman effect by in- 
troducing a gap in the energy dispersion, then calculate the electron and the 
current densities together with the electrostatic and electrochemical potentials 
as a function of the lateral coordinate. Moreover, the longitudinal resistance 
(R X x) is calculated for varying magnetic field. In our calculations we assume 
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a translational invariant sample in y direction, where an electron channel is 
formed in the interval —b<x<b, with b being the depletion length. The ge- 
ometry under consideration is the so-called floating gate structure [4]. Given 
the boundary conditions and external potentials one can solve the Poisson 
equation, starting from the zero temperature, zero magnetic field solution. 
The basic ingredients of the SCTFPA are the iterative solution of the electron 
density, 

n el (x) = J dED(E)f(E - V(x)), (1) 

where f{a) is the Fermi function, D(E) a relevant density of states, and the 
total potential energy, 

V(x) = V ext (x) + (2e 2 /K) ^ n el {x')K{x,x')dx' 

+g*a/j B B. (2) 

Here, the first term stands for the external potential composed of gates and 
donors. In the second (Hartree) term k is the dielectric constant of the material 
and K(x, x') is the solution of the Poisson equation preserving the boundary 
conditions V(—d) = V(d) = 0, where 2d is the sample width. The third term 
is the Zeeman energy with effective g factor (g* = 5.2 from experimental and 
theoretical findings [9] and the references therein), a = ±1/2 spin polarization 
and n B the effective Bohr magneton considering GaAs/AlGaAs heterostruc- 
tures. The well established wisdom is that, the exchange-correlation effects 
enhance the Zeeman splitting, however, direct numerical calculations of these 
potentials without taking into account screening lead to discrepancies in ex- 
plaining the experimental data. In this work we take the experimental value 
of the g factor without dealing with the difficulties that arise from explicit 
calculation of the effects of exchange and correlations [6,7]. 

The imposed current densities are described by a LOL, by making use of the 
translation invariance and equation of continuity, namely 

E y {x)=E° y , and j x (x)=0, (3) 

and 



J y (x)=E° y /p xx (x), E x {x) = P -^-El (4) 

Pxx \E) 

where p xx {x) and p xy (x) are the diagonal (longitudinal) and off-diagonal (Hall) 
components of (local) resistivity tensor, p{x) = [a^x)]" 1 , respectively. Hence, 
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Fig. 2. (a) The numerically calculated current distribution as a function of external 
magnetic field, lines present the equi-current contours, (b) R xx against magnetic 
field, obtained at default temperature. 

given the local conductivity tensor elements one can obtain the current density 
in x direction and thereby the global resistances. We choose a simple (gener- 
alized) description of the conductivity model, namely the Hall component is 
(e 2 /h)v(x), where u(x) is the local filling factor and the longitudinal 
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component is a xx = 



a xy [e+(l-is(x)) 2 /4], < u(x) < 1.5 

a xy [e + (2 - is(x)) 2 /4], 1.5 < u(x) < 2.5, (5) 

where e is a cut-off parameter defining the accuracy of the numerics (~2x 
10~ 6 ). The above scheme, enclosures the simplest set of assumptions to provide 
a qualitative understanding of the current distribution and the calculation of 
the global resistances and can be improved in many aspects, however, as a 
first attempt grasps most of the essential physics. 

Fig. 1 depicts the calculated local filling factors versus the lateral coordinate, 
when the z/(0) < 2 (black solid line) and 2 < z/(0) < 3 (red solid line). The 
ellipses encircle the region where the local filling factor is one (dotted) or two 
(dash-dotted). It is seen apparently that, the current (broken lines) is essen- 
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tially confined into a region where an IS exists. In the case of u(0) < 2, most of 
the current is flowing either inside the IS or in the very close vicinity, whereas 
for 2 < i/(0) < 3 beside the main two peaks near filling factor two, we observe 
that there exists two other peaks on both sides coinciding with the positions 
of filling factor one regions (highlighted by arrows). This essentially indicates 
that although most of the current is carried by the innermost IS, if an external 
current is imposed there is a possibility for electrons to flow from these very 
narrow local regions, due to a local minimum of the longitudinal conductivity. 
This finding coincides with the very recent model by I. Neder [10], considering 
a non-Gaussian noise to explain the visibility oscillations performed at the 
Mach-Zehnder interference experiments. In Fig. 2 we show the local current 
intensity ( function of magnetic field together with the longitudinal 

resistance (b). The evolution of the reminiscence IS with u(x) = 1 can also be 
observed in the left side of the current distribution, however, the longitudinal 
resistance vanishes within the numerical accuracy in this interval. If the ISs 
become large than the magnetic length or the Fermi wave length, the system 
exhibits finite R xx and the current is distributed all over the sample, otherwise 
R xx = 0. The bubble-like features seen at the centers of the sample is simply 
due to the conductivity model we consider. 

In summary, we have calculated the spatial distribution of the current carry- 
ing, spin-split ISs within a Hartree type approximation incorporating Zeeman 
splitting by an experimentally induced effective g factor. We have shown that, 
the formation of (spin resolved) incompressible strips dominate the current 
distribution and thereby defines the magnetic field interval where one would 
observe the IQHE. 
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